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Abstract 

An imaginary part of the false- vacuum energy density in a metastable system, 
i.e., the decay width due to quantum tunneling, might be reproduced by Borel 
resummation of vacuum bubble diagrams. We examine the convergence of 
this prescription in the Gaussian propagator model, in which the analytical 
expression of vacuum bubbles to the ninth order of loop expansion is available. 
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Usually, one uses the bounce calculus |T| to estimate the decay width of a metastable 
system due to quantum tunneling. This is consistent with the conventional picture that the 
quantum tunneling is a purely non-perturbative phenomenon. However, several methods 
have been proposed so far which might allow us to handle tunneling effects with using 
a conventional perturbation series. According to these approaches, a certain class of physical 
systems seems to have in the perturbation series the whole physical information, including 
higher order corrections to the leading bounce approximation. 

In |§ , Suzuki and the present author proposed such a perturbative method to reproduce 
an imaginary part of the false-vacuum energy density. The method is based on Borel re- 
summation m of vacuum bubble diagrams and is applicable to super-renormalizable field 
theoretical models. The validity of the method has been tested in one-dimensional (i.e., 
quantum mechanical) and two-dimensional O(N) symmetric A0 4 model 0. In the former, 
the imaginary part of the ground state energy was successfully reproduced in a high accu- 
racy. In the latter, unfortunately, a shortage of the perturbation series prevents us to draw 
a quantitative conclusion on the convergence of the method. 

In this note, we test the convergence of our proposal by applying it to the scalar 4 
theory with a Gaussian propagator |J. The Euclidean action is given by 



S[(j)}= I d D x 



(1) 



where A is D-dimensional Laplacian. This model has advantages for our purpose: Since the 
propagator in momentum space is 

G{p)=e-?\ (2) 

an evaluation of Feynman diagram is just Gaussian integrations and no UV divergences 
present. These facts allow the evaluation of Feynman diagrams to a relatively high order in 
an arbitrary dimensional space . This model also allows to carry out the analytical bounce 
calculus Q. Unfortunately, the kinematical meaning of this model is not clear due to the 
odd-looking "kinetic term," (actually, the classical equation of motion has an infinite order 



time derivative) and thus the meaning of "tunneling" or "decay" is not obvious. Nevertheless, 
we adopt this model as the testing ground to investigate the convergence of our prescription 
for Borel resummation of vacuum bubble diagrams. 

We shall consider a "metastable case" (g > 0) and study a sum of vacuum bubble 
diagrams: 

oo 

£(g) ~ £ c n g n . (3) 

n=0 

The actual numerical values of c n s are given in |J to n = 8 (to nine loops) from D = 
to D — 4 and we will use those values. 

In an analogy of the conventional field theory, we shall call 8(g) "the vacuum energy 
density." The imaginary part of the vacuum energy density is also evaluated in the leading 
order bounce approximation ||: 



where 



A n = - — , J - 1 



2 V 2 J ' V g 

M = 2 



2\ D/2+1"! £> 



3y j expE 2(D ), S2 ( g ) = | 2 ( ^_-;f ln(l-3^). (5) 

This formula tells us the leading large order behavior of the perturbation series (§) @ , or 
equivalently, the position and nature of the Borel singularity nearest to the origin of the 
Borel plane. 

In completely the same reasoning as we thus define the Borel-Leroy transform of the 
perturbation series c n : 

B( z ) = Y — ^ — (6) 

The original vacuum energy density is supposed to be reproduced by the Borel integral: 

g(D+i)/2 J 



1 POO 

^) = -77577w/ dze-l'2p+*>^B{z). (7) 
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The leading imaginary part (Q) for g <C 1 implies that there exists a singularity of B(z) 
at z — A Q . The argument of the gamma function in (0) was chosen to sufficiently weaken 
this singularity: The bounce singularity at z = Aq becomes a square root branch point, irre- 
spective of D. The integration contour of (|7|), therefore, must avoid the bounce singularity. 
This deformation develops the imaginary part and we select the contour along the upper 
side of the real axis. 

The singularity at z — Aq, however, limits the convergence radius of the series (H). 
Consequently, in order to carry out the integral (0) beyond z = A , we have to execute the 
analytic continuation of B(z). This continuation is possible only if we completely know the 
Borel function B(z) within the convergence circle (i.e., if all c n s are known). However, we 
have only a finite number of perturbation coefficients in actual. 

To avoid this difficulty, we use the conformal mapping technique. Introducing a new 
variable A instead of z, 



z = AA - 



A 



(8) 



(1 + A) 2 ' 

new integration contour is confined in the convergence circle on A plane. Indeed, the bounce 
singularity is transformed to A = 1 and new integration contour goes along inside the upper 
arc of the convergence circle. 

In this way, we have an imaginary part of the vacuum energy density expressed by the 
perturbation coefficients fl3|: 



lm£ (g) 

_ fAo 
V 9 



p 

(D+l)/2 



d9 exp 



with 



A 



g cos 2 (6/2) 



sin(0/2) ? . 

d k sm(k6), 



cos D+2 (#/2) 



(9) 



k=0 



d uk -n r(fc + n)(4A)) w 

" „=o l ] {k-n)\r(2n)T{n + {D + l)/2) { 



(10) 



where P denotes the order of approximation in this approach. 
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Now we present the result of eq. @. When choosing the presented range of the coupling 
constant g in Figures, we selected as the upper bound the point that the effective mass 
changes it's sign. To the one- loop order, we have 

EW U = 1 -5(4^ + 0(92) - 

Thus we choose as the upper bound of g, g = 2(47r) D / 2 (numerically, for D = 1, 2, 3 and 4, 
it is g = 7.09, 25.1, 89.1 and 315, respectively). This upper bound might be regarded as an 
order of magnitude estimation of the boundary of "tunneling region" ||. 

We plot eq. @ with P = 4, 5 and 8 as the function of g. The imaginary part is 
normalized by the result of the leading bounce approximation (£|). For D = (Fig. 1), 
all the perturbative coefficients in (§) can be analytically computed because the partition 
function is just an ordinary integral. Eq. (|9]) rapidly converges to a uniform curve as shown 
in Fig. 1 and our prescription certainly seems to give the correct imaginary part in D = 0. 

For D = 1 (Fig. 2), three curves almost coincide over a wide range of g. In particular, the 
impressive agreement of P = 5 and 8 indicates that the perturbation series to the fifth order 
is practically sufficient to reproduce the correct imaginary part. For D = 2 (Fig. 3), the range 
in which P = 5 and 8 agree becomes narrower and simultaneously the difference becomes 
larger. But still eq. @ displays a rapid convergence for all g and therefore P = 8 seems to 
indicate the correct value. The order of perturbation series required for a plausible evaluation 
of the imaginary part is higher than the cases of D = and 1. 

For D = 3 and 4 (Figs. 4 and 5), there is a wild oscillation in small g region. (For 
g = 0+ the correct value is 1 because the bounce result @ becomes exact.) In spite of this 
undesirable situation for the weak coupling, eq. @ seems converging in the strong coupling 
region. Therefore, in both cases of D = 3 and 4, the correct imaginary part in the strong 
coupling region seems to be smaller than the bounce result. 

In summary, our method applied to the Gaussian propagator model seems certainly 
working in lower dimensional cases, exhibiting a good convergence. The results in turn 
indicate the correct imaginary part of the vacuum energy density is much smaller than that 
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of the leading bounce approximation. On the other hand, the convergence in the weak 
coupling region becomes worse in higher dimensional cases. We guess this tendency holds 
in conventional field theories. 

I would like to thank Dr. H. Suzuki for numerous discussions and a careful reading of 
the manuscript. I am grateful to Prof. J. Zinn- Justin for kind correspondence. 
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FIGURES 

FIG. 1. The imaginary part of the "vacuum energy density" evaluated by our prescription (|9|) 
versus g: P = 4 (boxes), P = 5 (circles) and P = 8 (filled circles). The imaginary part is normalized 
by the result of the leading bounce approximation (Q). 

FIG. 2. Same as Fig. 1 but for D = 1. 

FIG. 3. Same as Fig. 1 but for D = 2. 

FIG. 4. Same as Fig. 1 but for D = 3. 

FIG. 5. Same as Fig. 1 but for D = 4. 



S 



1.2 



o 
H 



0) 



0.8- 



0.6- 



. 4 



0.2 




Fig.l 




Fig. 2 




Fig. 3 




Fig. 4 




Fig. 5 



